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Based on the Boltzmann-Enskog kinetic theory, we develop a hydrodynamic the- 
ory for the well known (reverse) Brazil nut segregation in a vibro-fluidized granular 
mixture. Under strong shaking conditions, the granular mixture behaves in some 
ways like a fluid and the kinetic theory constitutive models are appropriate to close 
the continuum balance equations for mass, momentum and granular energy. Using 
this analogy with standard fluid mechanics, we have recently suggested a novel mech- 
anism of segregation in granular mixtures based on a competition between buoyancy 
and geometric forces: the Archimedean buoyancy force, a pseudo-thermal buoyancy 
force due to the difference between the energies of two granular species, and two 
geometric forces, one compressive and the other-one tensile in nature, due to the 
size-difference. For a mixture of perfectly hard-particles with elastic collisions, the 
pseudo-thermal buoyancy force is zero but the intruder has to overcome the net 
compressive geometric force to rise. For this case, the geometric force competes 
with the standard Archimedean buoyancy force to yield a threshold density-ratio, 
R p i = pi/ps < 1, above which the lighter intruder sinks, thereby signalling the onset 
of the reverse buoyancy effect. For a mixture of dissipative particles, on the other 
hand, the non-zero pseudo-thermal buoyancy force gives rise to another threshold 
density-ratio, R P 2 ( > R p i), above which the intruder rises again. Focussing on 
the tracer limit of intruders in a dense binary mixture, we study the dynamics of 
an intruder in a vibrofluidized system, with the effect of the base-plate excitation 
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being taken into account through a 'mean-field' assumption. We find that the rise- 
time of the intruder could vary non-monotonically with the density-ratio. For a 
given size-ratio, there is a threshold density-ratio for the intruder at which it takes 
the maximum time to rise, and above(/below) which it rises faster, implying that 
the heavier (and larger) the intruder, the faster it ascends. The peak on the rise- 
time curve decreases in height and shifts to a lower density-ratio as we increase 
the pseudo-thermal buoyancy force. The rise (/sink) time diverges near the thresh- 
old density-ratio for reverse-segregation. Our theory offers a unified description for 
the (reverse) Brazil-nut segregation and the non-monotonic ascension dynamics of 
Brazil-nuts. 

Keywords: Granular mixture — Brazil-nut segregation — reverse buoyancy — non- 
monotonic rise velocity. 
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I. INTRODUCTION 



The phenomenon of segregation, in which a homogeneous mixture of particles of different 
species becomes spatially non-uniform by sorting themselves in terms of their size and/or 
mass, is ubiquitous in numerous chemical and pharmaceutical industries, dealing with the 
transport and handling of bulk granular mixtures 

HQ. 

For most industrial processes, 
it is required to maintain a homogeneous mixture during processing, with the segregation 



being the 'unwanted' phenomenon. Despite its unwanted consequences^ 



rowever , segregat ion 



id, 
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occurs spontaneously in driven granular mixtures 0, 4], |], [f], Q, [§], O, 

When a mixture of different particles is subjected to vertical shaking in a container, the 
larger particles (intruders) rise to the free surface. This is known as the Brazil-nut phe- 
nomenon (BNP), and is one of the most puzzling phenomena of granular materials research, 
still lacking a proper theoretical explanation. The early experimental investigations on the 



Brazil nut phenomenon 



flflfl 



were mainly concentrated on the effect of size of the 
intruder particle in vibrated systems. Besides experiments, computer simulations have been 

n u n n 

extensively used to study segregation [3|, 13, |8j, 19J. In fact, the current interest within the 
physics community on granular segregation was stimulated by the Monte Carlo simulations 
of Rosato and co-workers 3J. They explained BNP as a percolation effect since the smaller 
particles can easily percolate down to fill the void, created behind the larger particles due 
to external shaking, which, in turn, pushes the larger particle to the top. However, we need 
to point out that such percolation effects are likely to be active in a dense bed only under 
weak-shaking conditions. 



In many experimental setups convection is unavoidable 



ingredient to drive segregation [4, O, Q] • In the regime of convection-driven segregation |4j , 
it has been experimentally verified that the large particles rise with an upward-plume of the 
surrounding bed at the center of the container, but cannot sink to the bottom since they are 
unable to fit themselves in a narrow downward-plume near the side-walls. It has subsequently 
been verified, via careful MD simulations 8], that convection-driven segregation dominates 
(over percolation) in deep beds. In this scenario the intruder particles with different sizes 
rise at approximately the same rate 4|]. 

However, a regime where convection is negligible has also been found in some experiments 

, lu, LOf] - The most striking observation in this 



and seems to be a ke^ 



15, 



5] as well as in many idealized simulations 



3 



case is the dependence of the segregation-rate on the particle-size, where the ascending 
velocity of the intruder particle increases with its size It has been found that there is 

a threshold size-ratio for the intruder above (/below) which the intruders rise (/sink). 

A series of recent discoveries have shown that the Brazil nut phenomenon is more intricate 
than it might seem. The subject took a major s tep forward with the papers of Shinbrot 
and Muzzio 10] and Hong, Quinn and Luding In 1998 Shinbrot and Muzzio [T" 



discovered a reverse buoyancy effect: while a large heavy intruder rises to the free surface, 
an equally large but light intruder sinks to the bottom of the granular bed. Three years later 
Hong, Quinn and Luding ll| introduced a phenomenological theory for the reverse Brazil 
nut phenomenon (RBNP): a competition between percolation and condensation, depending 
on the size and the mobility of the intruder particles, could drive the intruders to sink to 
the bottom of the container and vice versa. The key idea of this theory is that a mono- 
disperse system of particles can be ffuidized above a critical temperature (T cr ), assuming 
that the system is shaken uniformly. Hence, if a binary mixture is vibrated such that the 
mixture temperature lies between the critical temperatures of two individual species, then 
one species will be in a ffuidized state and the other condenses at the bottom, leading to 
segregation. Immedia tely the subject attracted the attention of researchers who performed 
new experiments 
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As mentioned before, different mechanisms have been proposed to explain the segrega- 
tion phenomenon, for example, percolation j^], arching j^], convection 4, f], 8], inertia 10]. 
condensation H| and interstitial- fluid effect s [l2 | . Unfortunately while the observational 



evidence accumulates 



Q y, q, y, y, 0, y, y, h 



rfh 

unified theory for the dynamics of Brazil nuts [HI, 115 



1 
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relatively little work exists on a 



171 ]. Recently, Jenkins and Yoon 



17| developed a theory for the segregation of elastic particles using hydrodynamics equa- 



tions of binary mixtures. They investigated the upward-to-downward transition introduced 



m 



111 ] even though their prediction of the phase diagram for the BNP/RBNP-transition 



did not match with the simulation results of Hong et al. Interestingly, none of these early 
theories considered either the non-equipartition of granular energy (which is a generic fea- 
ture of granular mixtures) or the effect of external driving forces. Thus, it appears that a 
comprehensive theoretical description for the dynamics of Brazil nuts is still lacking. 

A hydrodynamic model to include the effect of the non-equipartition of granular energy 
was postulated in ref. [18J where a new mechanism for segregation due to buoyant forces 
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was introduced, drawing a direct analogy with the buoyancy forces in fluids. More recently, 
a minimal hydrodynamic model for segregation was outlined in ref. flo| . starting from 
the Boltzmann-Enskog-level continuum equations for a dense binary mixture of fluidized 
granular materials. The important effect of dissipation, which is responsible for the non- 
equipartition of the granular energy, was also incorporated. This latter work clearly shows 
how one could derive a time-evolution equation for the relative rise velocity of the intruder 
from first principles, thereby making the analogy with standard fluids transparent (which 
has a history of more than a century, with the seminal works of Stokes, Oseen, Boussinesq, 
etc). It was argued in ref. [3] that the BNP/RBNP segregation dynamics results from 
a competition between buoyancy and geometric forces. This analysis 1 19] appears to be 
compatible with the experimental observations reported by Breu et al. [l3| and also agrees 
with the molecular dynamics simulations reported by Hong et al. 111 ]. 



In view of the good qualitative a gree ment of our model 19| with the experimental ob- 



servations presented by Breu et al. 



13], it is worth attempting to construct a theory for 



segregation in granular materials using the kinetic theory of inelastic hard particles. We 
show that a granular hydrodynamic theory does indeed provide a good description for the 
dynamics of Brazil nut segregation. 

This paper is organized as follows. In Section [TT] we introduce the kinetic theory for 
inelastic hard particles, by summarizing the principal definitions of the model and presenting 
the balance equations for the mixture. The problem of the non-equipartition of kinetic energy 
is considered in III B 21 The main results of this paper are presented in Section IHIl where 
the detailed theory for segregation is outlined. The evolution equation that governs the 
segregation dynamics of intruders in the dense collisional regime is outlined in IIII Al We 
then comment on the origin of different 'gravitational' segregation forces derived from the 
model in Section IIII B 11 followed by a discussion on the origin of the unsteady forces that 
act on the intruders in Section IIII B 21 In Section IIVI we discuss the phase-diagram for the 
BNP/RBNP transition in three-dimensions (i.e. for spherical particles), thereby verifying 
that the competition between the buoyancy and geometric forces drives the segregation 
process in a fluidized mixture. In addition, we explain the reverse buoyancy effect of Shinbrot 
and Muzzio [l(3| in Section IIV Al In Section |V] we apply our theory to probe the dynamics 
of intruders in the tracer limit, thereby explainin g so me recent experimental observations 
on the non-monotonic rise-time (Mobius et al. 12j) of the intruder (with the density- 
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ratio) as well as the divergence of the rise(/sink) time near the BNP/RBNP transition. In 
Section IVl Al we provide a simple analytical explanation for the RBN effect, focusing on the 
Boltzmann (dilute) limit. In Section IVlBl we discuss the possible higher-order effects of the 
non-Maxwellian velocity distribution on segregation forces, along with suggestions on some 
experimental implications of the present work. The conclusions and the limitations of our 
theory are detailed in Section IVIII 



II. KINETIC THEORY OF GRANULAR MIXTURES 



We assume that granular matter can be described at the 'macroscopic' level by a set of 
continuous hydrodynamic equations as a fluid-mechanical medium. The balance equations 
and the constitutive relations can be derived from a kinetic theory description. For granular 
media the kinetic equations are modified to account for the inelastic nature of the collisions 
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22] . 



between particles |14j]. These equations have been extended to binary mixtures |20l . 
The validity of the hydrodynamic approach even in the dense limit has recently been justified 
via a comparison of the theory with various experiments 231 ] . 

Theoretical and numerical studies for binary granular mixtures have shown that the two 
components have different kinetic (fluctuation) energies 



20 



22 



24 



25 



261 ] which has also 



been confirmed in vibrofluidized experiments [271 ] . Recently, Alam and Luding [24], [25] 
have demonstrated that a proper constitutive model for granular mixtures must incorporate 
the effect of the non-equipartition of granular energy. Therefore, in the present model we 
consider this breakdown of equipartition on the kinetic energy and this is a fundamental 
difference with previous related studies [ll|, Q, 12 1- 

The following subsection intends mainly to establish some notational conventions used 
throughout the paper and recall some basic definitions on granular hydrodynamics. Here 



we follow the lines of refs. 



20 



24 



2l| in the introduction of the kinetic theory model and refs. 



251 ] in the discussion of the breakdown of the equipartition of the kinetic energy. 



A. Definitions 



As a mechanical model for a granular fluid we consider a binary mixture of slightly 
inelastic, smooth particles (disks/spheres) with radii r\j (i = l,s, where I stands for large 
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particles and s for small), mass rrii in two or three dimensions (d = 2, 3). The coefficient of 
restitution for collisions between particles is denoted by e^, with < 1 and = e^. 

The average macroscopic quantities are calculated by taking appropriate moments of the 
corresponding microscopic (particle-level) property in terms of the single particle velocity 
distribution function fi(c,r,t) for each species. By definition fi(c, r, t)dcdr is the total 
number of particles which, at time t, have velocities in the interval dc centered at c and 
positions lying within a volume element dr centered at r, i.e., 

J fi(c,T,t)dcdT = N t . 

If the particles are uniformly distributed in space, so that fa is independent of r, then the 
number density rii of species % is 

rii = J /i(c,r,£)dc, 

and the total number density is n = rii + n s . The species mass density Qi is defined by the 
product of rii and rrii, and the total mixture density is 

Q = Qi + Qs = Pi4>i + Ps4>s, 

where pi is the material density of species i and <pi is the ci-dimensional volume fraction for 
species i: 

= ~d~ UiTi ' 

where Qd is the surface area of a <i-dimensional unit sphere. 

The mean value of any quantity ipi = ip(c) of a particle species i, is 

<^(c)> = - f ^(c)/,(c)dc. 

The mean velocity of species i is Uj = (cj). The mass average velocity u of the mixture is 
defined by 

u = -{giui + g s u s ). 
Q 

The peculiar (fluctuation) velocity of species i is Cj = Cj — u, and the diffusion velocity is 
Vj = (Cj) = Uj — u. The species "granular temperature" is defined proportional to the mean 
kinetic energy of species i 

Ti = ~^ n i (C« • Cj) , 
7 



and the mixture temperature is 

T = -(niTi + n s T s ). 
n 

Let us remark that this generalized notion of temperature is introduced for a theoretical 
convenience to take advantage of a thermodynamical analogy for granular materials, and 
thereby postulating a higher-order field variable. Even though the definition of thermody- 
namic variables for non-equilibrium states is straightforward theoretically, the thermody- 
namics of non-equilibrium states has always been a matter of debate which we will not touch 
upon here. 



B. Mixture Granular Hydrodynamics 

The evolution of the granular system is governed by the well known balance equations 
for the mixture density, momentum and energy: 

g = -gV-u, (1) 

gu = -V -P + ^mFi, (2) 

i=l,s 

= TV ■ j - V • q - P : Vu + £ j< • F 4 - V, (3) 

i=l,s 

where the overdot indicates the convective derivative: d t (-) +u- V(-). Here, P is the mixture 
stress tensor, Fj is the external force acting on the particle, q is the mixture energy flux, j 
is the diffusive mass-flux and T> is the total inelastic dissipation rate. These equations are 



rigorous consequences of the Enskog-Boltzmann kinetic equation 28j . extended to inelastic 



particles [20|, |2l|], and must be supplemented with respective constitutive equations for P, 



q, j and T>. For = 1, the collisional dissipation rate vanishes [T> = 0), and consequent 



28) 



y 



we recover the standard energy balance equation for a mixture of elastic hard spheres 

1. Species Balance Equations and Fluxes 
The species momentum balance equation is: 

I Ti ■ 

dtiriiUi) + V • (riiUiUi) = V • Pi + —F { + T h (4) 

rrii rrij 



where I\ is the momentum source which arises due to the interaction between unlike particles 
and Y2i=i s = 0- 111 the limit of small spatial inhomogeneities (i.e. first order in the 
gradients of the mean fields) the species stress tensor Pj has, at the Navier-Stokes level, the 
standard Newtonian form 

Pi = p i I + // i (Vu + Vu T ), (5) 

where pi is the partial pressure of species i, \ii is the viscosity of of species % and I the uni t 
tensor. The equation of state for the partial pressure of species % can be written as 
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211: 



Pi 



^i^i^ii 



where is the "compressibility" factor of species i, 

Z, = 1 



j=l,s 



(6) 



(?) 



and 



1 



K a = o&XijCt + Ri. 



Here Xij is t ne contact value of the radial distribution function and Rij 
ratio. The radial distribution functions o 



(8) 

ri/rj the size 



granular systems are often app roximated by their 
elastic counterpart (see, for example . |20l . \21\\ and for related issues [29]). In this paper we 
use the following functions for disks 3( 

Xij 

and for spheres: 



1 
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nRu + <j) s Ri, 



1 



8(l + ^)(l-0) 2 



(9) 



1 



Xij 



1 - 



+ 



nRu + cj) s Ri 



3 + 



>iRu + cj) s Ri 



(10) 



An explicit expression for the momentum source term can be written as [20] ]: 

mk-mi" 



riiK ik T 



4 



rriik 

o \ 1/2 

zmiirik \ 



V(lnT) + V 



Inl^l 



(11) 



r ik \nm ik T / 

where = m, + and = r, + r^, with i ^ k. In deriving the above expression, 
it has been assumed that the single particle velocity distribution function of species i is a 
Maxwellian at its own granular energy Tf. 



This represents the zeroth-order approximation for the distribution function which has re- 
cently been verified in molecular dynamics simulations of vibrofluidized binary granular 
mixtures 3jJ. Higher-order corrections to the distribution function would appear beyond 



the Euler- level description that we neglect in the present work (see discussion in VLB). 
The energy balance equation of species % is 

dtintTi) + V • {uiUiTi) + V ■ Qi + Pi : Vu 

-mFi ■ v, - ^v, • (V ■ Pi - nF) = T>i, 
Q 



and the mixture energy flux is defined as 



q = E HiTi + 



i=l, 



where Qj is the energy flux of species i: 



(13) 



(14) 



(15) 



and Ki is the analog of thermal conductivity of species %. The rate of kinetic energy dissipa- 
tion of species i, Vi, is [26] 

(16) 



V 



j=l,s 



Ti Ta 



m,- m* 



T- - T 

rrii 



Ti 

— + — 

VTLi rrii 



1/2 



3 / " "V 

where Mji = mj/mij. The total rate of kinetic energy dissipation is simply T> = T>i + T> s . 
We can split eqn. ffTB"]) into two terms: 



1. The inter-species collisional dissipation rate, V(, 



2. The exchange collisional dissipation rate, Vf 



2V2d 



Xl kR? k (l + RkiT^M^M^l + e lk ) 



„d+l 1/2 

r i m/ 



T 



1 + 



mT k \ 

m k TiJ 



(18) 



with k 7^ i. 
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Note that J2i=i s = 0; an d the exchange term Vf is a consequence of the non- 
equipartition assumption. On the other hand, with the equipartition assumption (T} = 
T S = T), Vf = and Eq. Qnj reduces to: 



2. N on- equipartition of Granular Energy 

It is worthwhile now to comment on the breakdown of the equipartition partition principle 
for a granular mixture. Physically the lack of energy equipartition is determined by the 
different dissipation rates ifPTj) and ( 1181) . An obvious question is, "How each species will 
attain a different granular temperature?" The answer to this question depends on the 
balance between the external power injected into the system and the energy dissipation 



rates for the two species [22|, |24], |25|, |26|], for example, in a vibrofluidized bed. 

Before formulating a theory for segregation, it is important to have an estimate of the 
granular energy ratio T\jT s in a fluidized granular mixture. An adequate description of the 
vibrofluidized granular system requires a careful analysis of the steady version of the species 
energy balance equations as described above. To formulate the associated boundary-value 
problem, however, we need to impose appropriate boundary conditions at the base-plate 
and the free surface of the vibrofluidized mixture. The problem of boundary conditions to 
be satisfied for the granular hydrodynamic equations of a vibratory mono-disperse system 



has been studied [32|], but such studies for binary and polydisperse mixtures do not exist. 
Hence the complete solution of the boundary-value problem for a vibrofluidized mixture is 
left to a future effort. 

In the results presented below, we have rather used the granular energy ratio of Barrat 
and Trizac [26j which was derived for a randomly heated granular gas (see Section I V) . 



Their expression has subsequently been verified in the molecular dynamics simulations [31] 
on vibrated granular mixtures. Paolotti et al. [31] showed that the granular energy ra- 
tio remains constant in the bulk of the mixture, with the variations of T;/T s with height 
being concentrated in two boundary layers near the base plate and the free surface. The 

n 

experiments [27| also showed the existence of a constant Ti/T s in the bulk. 
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III. THEORY OF SEGREGATION 



To illustrate how the segregation dynamics can be described by a hydrodynamic approach, 
we make the following approximations: 

• We assume that the granular medium is in a fluidized state. The fluidized state can 
be realized when the granular material is vibrated strongly in the vertical direction, 
typically by a harmonic excitation z p (t) = A sin (ut), with the amplitude A and the 
frequency uj = 2nf. The normalized acceleration parameter (r = Au 2 /g) satisfies 



• We assume that the system is in a regime where the bulk-convection and the air-drag 
can be neglected. 

• We neglect the viscous stresses (as well as any stress-anisotropy) for the case where 
there is no overall mean flow in the system. For example, in a typical experimental 
realization of a vertically vibrated bed the mean velocity is zero (u = 0). 

• We impose horizontal periodic boundary conditions to make the equations analytically 
tractable. Therefore, the hydrodynamic fields vary only along the vertical direction 



Here we outline the derivation of the evolution equation for the Brazil nuts. With the 
above assumptions, the equations for the evolution of the granular mixture, ([I])-®, (TjJ 
and (TL3I) . simplify considerably. Our starting point here is the inviscid species momentum 
balance equations (jlj): 



r > 1. 



(d/dz{-) ^ 0, but d/dx{-) = and d/dy{-) = 0). 



A. Segregation Dynamics: Evolution Equation 



d t ui = d z pi - g H F h 

01 01 



(19) 



d t u s = d z p s - g H r s . 

0s 0s 



(20) 



Subtracting Eq. ( 1201) from Eq. ( |T9l) we obtain 
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where u\ = u\ — u s is the relative velocity of the larger particles. Now we need expressions 
for d z pi and d z p s . 

From the equation of state (jSJ), we define the following weighted-ratio of two partial 



pressures 



0: 



9(0;, fa/fa, n/rj, Ti/Tj 



Ps/jh 

Pi/ni 



Pi n s 



which depends on the species volume fractions, the volume fraction ratio, the size-ratio and 
the temperature ratio. The partial derivative of pi is calculated from 



d z pi 



d (p s ni\ p s d fni 



dz \ O n. 



O dz V n 



+ 



1 dp s _ p s dQ \ 

e dz e 2 dz J 



Substituting this in Eq. (l2Tj) and using the expression for d z p s in Eq. (f2"0l) (and after some 
algebraic manipulations), we arrive at the following evolution equation for u\: 



Qid t Ui 



m. 



-Pid z 







\Z S 


S 


In ( 


ni 




n s Q 



+ 



1 + 



+ ni 



ZiTt 



d t u s . 



(22) 



In the following we focus on the tracer limit of this equation for a mixture where the 
number density of intruder (larger) particles is much smaller that of the smaller particles, 
i.e. n\ << n s . In this tracer limit (n; << n s ), the intruders are assumed to stay far away 
from each other (this, of course, implies that larger and smaller particles are homogeneously 
mixed initially) and they do not influence each-other's motion; in fact, as explained by 



Lopez de Haro and Cohen [331 ]. except for the mutual and thermal diffusion coefficients, the 
transport coefficients of binary mixtures with one tracer component can be obtained from 
the general expressions of corresponding transport coefficients by taking the limit n\jn s — > 0. 

For analytical simplicity, we make two assumptions: (1) the global temperature of the bed 
T does not vary with height (see, for example, the experiments in 34(); (2) the temperature 
ratio Ti/T s remains constant in the bulk which has also been verified in experiments 27| as 



well as in simulations 



3l| of vibrated granular mixtures. In the dense collisional regime of 



a mixture with n\ << n s and d z T = 0, the term associated with the momentum exchange 
in Eq. ( 1221) can be approximated by: 



1 + 



ni 

n s <3> 



T 



lnl^ 



4 / 2mim s 



ris \vrm Zs T 



1/2 
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Furthermore, in this limit, it can be shown that Zi/Z s ~ 4> 2 /2 and hence d z \n(ZiTi/Z s T s ) ~ 
d z ln<j) 2 = 2A, where A is the decay rate of the volume fraction with height: ln</> = 5 + Xz. 
An approximate value for A can be estimated from the vibrofluidized experiments 341 ] : 
Xr s < —0.05 (see Fig. 4c in 34j) which is assumed to be valid in the dense regime. Now, 
the third term on the right-hand side of Eq. (j22p and the first term in the above expression 
can be combined to yield 



Pi 



3 2 



In 



Pid z 



In 



ni 

n.,9 



Pi 



13, 



In 



4n 



+ Pid z 



In 



Z a T s 



Hence, both these terms are of order O(Xpi). With these assumptions and retaining terms 
of same order in Eq. (122]) . the time-evolution equation for the relative velocity of intruder 
particles u\ is described by the following equation: 



mi 



dt 



ZiT;, 
Z s T s 



- mi 



4:Ki s T ( 2mim 



ris 



nmi s T 



1/2 



+ 



m. 



ZhTj, 
Z R To 



mi 



du s 
dt 



(23) 



where all expressions are evaluated in the limit n\/n s — > 0. 

The above equation (l23jl contains all the information necessary to describe the segregation 
dynamics of intruders in the tracer limit. The first term on the right hand side is the net 
gravitational force acting on the intruder. Note that the second term has a form similar 
to the Stokes' drag force, which always acts opposite to the intruder's movement. The 
functional form of this term can be justified by recalling the experimental observations of 



Zik et al. 



35], who found a linear dependence of the drag force on the velocity of a sphere 



moving in a vibrofluidized granular medium. This could be interpreted as a characteristic 
of the fluidized state of the granular media. The last term represents a weighted coupling 
with the inertia of the smaller particles. In the following subsection, we discuss the origin 
of these forces in detail. 



B. Segregation Forces 

Let us now analyse different forces that are acting on the intruder as it rises/sinks through 
the granular bed. 
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1. Gravitational forces: buoyancy and geometric forces 



We can decompose the net gravitational force in Eq. (|23|) on an intruder in the following 
manner: 



(i>, - i>i)Vi + /». [y { jy- "'- (' ~ "rj ' " ls \ Y 1 



(24) 



where is the volume of a particle of species i. This net gravitational force is composed of 
the following buoyant and geometrical forces: 

• An Archimedean buoyancy force due to the weight of the displaced volume of the 
intruder (Vj): 

F b A := V^ps - fn)g. (25) 

• An analogue of the thermal buoyancy force due to the difference between the two 
granular temperatures: 

if :=/3(T,-T s ), (26) 

where j3 = m s T~ 1 {Zi/ Z s ) is the effective pseudo-thermal expansion coefficient. We 
call this pseudo-thermal buoyancy force. Note that this force vanishes identically for 
a mixture of particles with elastic collisions (e^ = 1). 

• Due to the size-disparity between the intruder and the smaller particles, the intruder 
has to overcome a compressive volumetric strain, e s * := (Vi/V s — 1). This results in a 
static compressive force of the form: 

J# := -m s ge s v \ (27) 

This force is always negative since the intruder has to rise against gravity. 

• A dynamic tensile force from the pressure difference due to the interaction between 
the intruder and the smaller particles: 

:= m s gef\ (28) 

where e^ yn := (Zi/Z s — 1) > can be associated with a weighted volumetric strain, 
tensile in nature. 
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The last two forces (1271) and (125]) are not related to standard buoyancy arguments. Thus, 
purely geometric effects due to the size- disparity contribute two new types of segregation 
forces: 

F ge = F°l + F*» = ~m s (e? - g. (29) 

Overall, the collisional interactions help to reduce the net compressive force that the intruder 
has to overcome. 

It is interesting to find out whether we could get back the standard Archimedes law from 
Eq. (1241) . This corresponds to the case where a large particle is immersed in a sea of small 
particles with r\ » r s . For this limit, it follows that 

F£ n - m.(yi/V. - 1) = -Ff e (30) 

and the net geometric force is 

F ge = 0. (31) 

Thus, the net gravitational force on a particle falling/rising in an otherwise quiscent fluid 
(at the same temperature) boils down to the Archimedean buoyancy force: 

F = F$ = g(p s - Pl )V t . (32) 



19] 



2. 



Unsteady forces: the added mass effect 



We have noted in the previous section that the inertia of the smaller particles is directly 
coupled with the motion of the intruder in the evolution equation (viz. eq. [23]) . For this 
unsteady (inertial) force t F also, we follow our earlier decomposition: 



du s 
dt 



(p s - PiWi +m s [ -J- - 1 | -J- + m s ( 1 - ^ | +m s ( -4- - 1 



X - J z s ' V" vj ' " \z s 1 (oo) 

It is evident now that the net unsteady term has contributions from the standard added 
mass force along with two new forces as described below. 

• As in the unsteady-motion of a particle in a fluid, the intruder in a granular mixture 
has to rise along with its surrounding smaller particles. The 'added' inertia of the 
displaced smaller particles that are being carried by the intruder gives rise to an 

on the intruder: 

due 



effective added mass force 



Vi(p s -pi) 



dt 



(34) 
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As expected, this force vanishes if the material density of the intruder is the same as 
that of smaller particles. 

• A thermal analog to the added mass force, due to the difference between the two 
granular temperatures, is given by: 

'^--KHt- (35) 

This force vanishes identically for T\ = T s that holds if the particle collisions are 
perfectly elastic. Hence this represents a new force for the granular system. 

• Lastly, we have exact analogues of the two geometric forces as described in the previous 
section: 

While the former is the analog of static 'compressive' geometric force, the latter is the 
analog of the dynamic 'tensile' geometric force as discussed in the previous section. 
Both these forces vanish if the intruder is of the same size as the bed-materials. 

Thus, we have shown that the last term in the evolution equation (1231) can be represented 
by an weighted added-mass force. 



IV. PHASE DIAGRAM FOR BNP/RBNP TRANSITION AND THE REVERSE 

BUOYANCY EFFECT 



To proceed in the simplest possible way, we first consider the steady-state solution of 
equation (1231) - in this case the added-mass forces do not influence the onset of segregation. 
Neglecting transient effects, the steady relative velocity of the intruder can be obtained from 

rug ( nmis • 
AK is \2mim s T j 



u, 



m* 



Zl_T\ 
Z K To 



mi 



(38) 



Setting this relative velocity to zero, we obtain the criterion for the transition from BNP to 
RBNP: 



17 



which agrees with the expression of Jenkins and Yoon [17| for the case of equal granular 
energies (T; = T s ), i.e. Zi/Z s = mi/m s . We have already noted in our previous paper 
[19I that the non-equipartition of granular energy (TJ 7^ T s ) that arises from the dissipative 
nature of particle collisions must be incorporated into the theory to correctly describe many 
experimental findings. 

As mentioned before, the energy ratio, Rt = Ti/T s , is calculated from the model of Barrat 



and Trizac 



26|]: 



C X R% 2 + C 2 ( 1 + —R T ) + C 3 ( 1 + —R T ) (R T - 1) + C 4 = 0, (40) 
mi I V m i 



3/2 / \ 1/2 



where 



/ \ 3//2 

C 3 = 2v / 2(l + e Zs )(l + i? s/ ) d - 1 M s/ (0 s M si + /j Rf / M is )x^, 
Ca = _2 d - 1 (l-eL)0 s ^r 1 X,s. 

To draw the phase-diagram in the (mi/m s , r;/r s )-plane, we need to solve to the segregation 
criterion ( |39l) in conjunction with the expression for T\jT s ( ]40l . This leads to a quadratic 
polynomial for the mass-ratio mi/m s , resulting in multivaluedness for mi/m s below a critical 
value of the size-ratio ri/r s as we discuss in the following. We note that another choice of 
Ti/T s , as suggested by one of the referees, does not change the qualitative nature of our 
results. 

In Fig. [T] we plot the phase diagram for the BNP/RBNP transition for a mixture of 
spheres in the tracer limit (4>i/<f) s = 1CT 8 ) at a total volume fraction of = 0.5. Note that 
this volume fraction is well below the value that corresponds to the perfect cubic-packing 
(0 = 7r/6 ~ 0.52), implying that the mixture is in the 'liquid' regime. Each solid curve in Fig. 
[U for a specified restitution coefficient, demarcates the zones of BNP and RBNP transition. 
We observe that the qualitative nature of the phase-diagram changes even if the particles are 
slightly inelastic (e = 0.99); for example, the mass-ratio (mi/m s ) is a multi-valued function 
of the size-ratio for e 7^ 1, in contrast to the perfectly elastic case (denoted by the dashed 
line). The effect of dissipation is to introduce a threshold size-ratio above which there is 
no RBNP. Moreover, the zone of RBNP shrinks dramatically when the particles are more 
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dissipative; for this parameter combination, there is no reverse segregation for moderately 
dissipative particles e < 0.8. We shall come back to discuss this point in the next subsection 
in connection with the reverse-buoyancy effect lo| . 

As discussed in our previous paper [19J, we need to look at various segregation forces to 
understand the driving mechanism for BNP/RBNP transition. First, we focus our attention 
to a mixture of equal density particles (pi = p s ); this case is easily amenable to experiments 
by using the intruder and the smaller particles of the same material. 

For an equal density mixture, the variations of different segregation forces with the size- 
ratio are plotted in Fig. El The coefficient of restitution is set to 0.95. Since pi = p s , the 
Archimedean buoyancy force, F b A oc (pi — p s ), is identically zero as shown by the dotted 
horizontal line in Fig. [2J The upper inset in Fig. [2] shows the variations of two geometric 
forces, and the net geometric force, F ge , is negative as shown by the dot-dash line. How- 
ever, the pseudo-thermal buoyancy force, F b oc (T} — T s ), is positive and increases with 
increasing size-ratio. Thus, the net gravitational force, F = Fb + F ge = F b T + F ge , can be 
positive/negative when F b T greater/less than F ge , respectively, with the equality being the 
BNP/RBNP transition point. The solid line in Fig. [2] shows the variation of the net grav- 
itational force F that changes sign at a size-ratio of Ti/r s ~ 2.6 above which the intruder 
rises (BNP) and below which it sinks (RBNP). Clearly, this threshold size-ratio is decided 
by a competition between the buoyancy forces and the geometric forces, leading to the onset 
of BNP/RBNP transition [lij. Note that for a mixture of particles with elastic collisions 
(e = 1), the pseudo-thermal buoyancy force is F b T = and hence F = F ge < 0; the intruder 
in such a mixture will, therefore, sink for the parameter combinations of Figs. Q] and [2j 

Now we focus on a mixture with the mass-ratio, m;/m s , between the intruder and the 
bed-particles being fixed; this can be realized in experiments by fixing the bed materials 
(glass beads or steel balls, etc.) and subsequently by varying the size and the density of 
the intruder. The variation of segregation forces with the size-ratio is shown in Fig. [3] for 
mi/m s = 10; other parameters are as in Fig. [TJ Both the buoyancy and geometric forces 
are negative at Ri s = 1 as seen in Fig. [3b,. While the net buoyancy force increases with the 
size-ratio, the net geometric force decreases in the same limit, and again the competition 
between these two forces decides the onset of BNP /RBNP. We observe in Fig. [3f? that the 
pseudo-thermal buoyancy force F b T remains positive only upto a size-ratio of Ri s < 11.8 and 
becomes negative thereafter (this is a consequence of the energy-ratio T\jT s being less than 
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1 for Rig > 11.8, see inset in Fig. (3j>). For this mixture, however, the Archimedean buoyancy 
force overwhelms its pseudo-thermal counterpart (i.e. F^ » F^) beyond a moderate size- 
ratio of Rig > 3. Hence, it is the competition between the Archimedean buoyancy force and 
the geometric forces that mainly determines the threshold size-ratio for the BNP/RBNP 
transition (Ri s pa 2.5) for this case. 



A. Reverse Buoyancy and Beyond 

In 1998 Shinbrot and Muzzio 



10] discovered an interesting effect: even though large 
heavy intruders can rise to the top in a vibroffuidized mixture of smaller particles, relatively 
lighter intruders of the same size can sink to the bottom. This is in contradiction to com- 
mon expectation and has been appropriately dubbed the reverse buoyancy effect. (Strictly 
speaking, however, the rising-phenomenon of 'heavier' intruders is also a reverse-buoyancy 
effect.) In the following we refer to the sinking-phenomenon of 'lighter' intruders {pi/p s < 1) 
as the reverse-buoyancy effect as in Shinbrot and Muzzio. 

To explore whether our model is able to explain this effect or not, we have fixed the 
size-ratio at Ri s = 2.5, and changed the mass-ratio by varying the density-ratio; again this 
can be realized in experiments by varying the material density of either the intruder or 
the smaller particles. Figure shows the onset of the reverse buoyancy effect where we 
have plotted the variations of different segregation forces with the mass-ratio, mi/m s , with 
other parameter values as in Fig. [2j To contrast the mass-ratio and density-ratio effects 
simultaneously, the same figure is redrawn in Fig. Hfr, with the abscissa now representing 
the density-ratio, pi/p s - For a fixed size-ratio, the net geometric force remains constant as 
shown by the dot-dash line in Fig. HI The upper inset in Fig. H] shows the Archimedean and 
pseudo-thermal contributions to the total buoyancy force; while F 6 T increases with the mass- 
ratio, F^ decreases in the same limit, and the total buoyancy force varies non-monotonically 
with both mass- and density-ratios. 

Focussing on the mass-ratio of mi/m s = 30 (the density-ratio is pi/p s ~ 1.98) in Fig. 
IU we note that the net gravitational force is positive, i.e. the heavier intruder will rise to 
the top (which is the BNP). Looking at the upper inset in Fig. [4h, we find that it is the 
pseudo-thermal buoyancy force that drives the intruder to the top (at such large density- 
ratios). Now if we decrease the mass of the intruder (by decreasing its material density p{) 
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to mi/m s m 25 (that corresponds to a density-ratio of pi/p s ~ 1-6), the net gravitational 
force becomes negative. Hence this relatively lighter intruder (pi/p s < 1.6) will now sink to 
the bottom (i.e. which is the RBNP). When we lower the density-ratio even below unity, 
the lighter intruder sinks to the bottom as observed in Fig. [4J. This is nothing but the 
reverse-buoyancy effect of Shinbrot and Muzzio [h]]: although the heavier intruder can rise 
to the top, equally large but a relatively lighter intruder can sink. 

It is interesting to observe in Fig. H] that there is a window of mass- and density-ratios 
(10 < mi/m s < 25 and 0.64 < pi/p s < 1.6) for which the net gravitational force remains 
negative and hence the intruder will sink. For pi/p s < 0.64, however, the Archimedean 
buoyancy force prevails over other forces and the lighter intruder rises to the top (as observed 
in the upper inset of Fig. 0J). This implies that if the material density of the intruder is 
much less than that of the bed particles, the intruder will eventually show the standard 
buoyancy effect (which is again the BNP). 

To understand the origin of reverse buoyancy, we consider particles with perfectly elastic 
collisions (e = 1), for which T\ = T s and _F fe T = 0; hence the net gravitational force, F = 
+ F ge , will decrease monotonically with increasing size-ratio as shown by the red-line 
in the lower inset of Fig. In such a mixture of elastic particles, an intruder with 

Rpi — Pi I ' Ps > 0.44 will sink (i.e. the reverse buoyancy), and this threshold density-ratio is 
less than unity since the lighter intruder has to overcome the net compressive geometric force. 
Thus, the onset of reverse buoyancy results from a competition between the Archimedean 
buoyancy force and the net compressive geometric force. 

The effect of the pseudo-thermal buoyancy force (for dissipative particles) is simply to 
increase this threshold density-ratio to R p i ~ 0.64 for parameter combinations of Fig. HI 
and create another threshold density-ratio at R P 2 ~ 1.6, above unity, beyond which the 
heavier intruder rises to the top (i.e. BNP). Thus, the lower threshold density-ratio R p \ is 
created by the geometric forces and the higher threshold density-ratio R P 2 is created by the 
pseudo-thermal buoyancy force. 

It is interesting to point out that the threshold density-ratio R P 2 at which the RBNP 
occurs depends crucially on the overall mean volume fraction (that can be related with the 
shaking strength of the vibrator, see Section IIVBI) of the fluidized bed. For example, by 
reducing the mean volume fraction, R p2 can be pushed below unity. 

Thus, our model predicts that the intruder will show the BNP for pi/p s » 1 and 
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pi/ p s « 1- While the BNP in the former limit is driven by the pseudo-thermal buoy- 
ancy force, the BNP in the latter limit is driven by the Archimedean buoyancy force. For 
intermideate values of density-ratios, the intruder will show the RBNP (and the reverse 
buoyancy) which is driven by a competition between the buoyancy and geometric forces. 

Now we make a few comments on the experimental realizations of reverse buoyancy. The 
range of size-ratios for which our model predicts the reverse buoyancy effect is of order 
Ris = 0(5) for the coefficient of restitution e < 0.99; for larger values of Ri s our model 
predictions show the standard Brazil-nut segregation. However, this critical size-ratio can 
be pushed to a much higher value by considering higher values of e; this is equivalent to 
reducing the magnitude of the pseudo-thermal buoyancy (F^), since F b T vanishes as the 
energy ratio Ti/T s approaches its equipartition value for perfectly elastic system (e = 1). 



More recently, Yan et al. 



37| | have also observed the reverse buoyancy effect for a range of 



size-ratios varying between 25 to 36; the diameter of their bed-particles was less than 250pm. 
Interestingly, the critical density at which this transition occurred was R p i = pi/p s ~ 0.7. 
But the interstitial air-pressure has played a crucial role in their experiments since they 
could not observe RBNP when the experiments were performed at a reduced air-pressure (0.1 
atm). Thus, the missing link could be provided by the air-drag that we have neglected in 
our analysis. 

Our predictions of a transition to BNP for very light intruders (at R p « R p i) has been 
recently observed by Huerta and Ruiz-Suarez 39] . Their experiments at high frequencies 
(/ = 50 Hz and T = 3) would closely mimic most of our assumptions since bulk-convection 
was negligible in their experiments. More careful experiments at vacuum are needed to map 
out the whole phase-diagram. 



B. Effect of Shaking Strength 

Here we consider the effect of the shaking strength of vibration on the BNP/RBNP 
transition. In typical vibrofluidized-bed experiments, the mixture is vibrated via a sinusoidal 
excitation of the base-plate: 

Zp{t) = A sin(27r ft), (41) 
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where A is the amplitude of vibration and / is its frequency. The shaking strength of 
vibration can be measured via the following non-dimensional number 

r = — i — J -L. 42 
9 

Ideally, the fluidized regime corresponds to both large-amplitude and high-frequency vibra- 
tions which can be related to the case T >> 1. Note that this condition, T >> 1, can be 
realized either by increasing the shaking amplitude at fixed / or by increasing the shaking 
frequency at fixed A. 

In the fluidized-regime, the bed expands with increasing shaking strength T, leading to a 
decrease in the overall volume fraction of the mixture <fi. Thus, the effect of V on the phase- 
diagram for BNP/RBNP can be tied to the effect of varying the mixture volume fraction 0. 
This is shown in Fig. [5] where we have plotted three curves for different volume fractions 
(0 = 0.1,0.3,0.5) with <f)i/(f) s = 10~ 8 . We observe that the range of size- and mass-ratios, 
for which the RBNP exists, increases with decreasing 0. This implies that the possibility of 
RBNP will increase with increasing shaking strength T. The experimental findings of Breu 
et al. [l3| showed similar trends with T, which our model is able to mimic. 

To understand the effect of mean volume fraction on the BNP/RBNP transition, we 
analyse different segregation forces for a given size-ratio and mass-ratio (say, at r\jr s = 4 
and rrii/m s = 40 in Fig. [5]). For this, both the Archimedean buoyancy force and the 
static geometric force F^ e do not vary with 0. The dynamic geometric force F^ 71 , however, 
decreases with decreasing 0. In this limit, the collision frequency decreases that decreases 
the collisional pressure, and hence the intruder will become relatively less mobile. For this 
case, even though the granular energy ratio do not vary appreciably with 0, the pseudo- 
thermal buoyancy force F^ will decrease with as a result of a decrease in the value of 
Zi/Z s . Both these effects result in a lower value of the 'upward' force on the intruder and 
hence the intruder will show the RBNP with decreasing mean volume fraction (i.e. with 
increasing shaking strength). For these parameter conditions, the mean density at which 
the BNP/RBNP transition occurs is about 0.284. 
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C. Effect of Intruder Volume Fraction 



So far we have shown results for a binary mixture in the tracer limit of intruders (i.e. 
( t ) i/ ( t ) s — 10~ 8 ). Recall from Fig. [T]that the zone of RBNP shrinks to zero as we make the 
particles more and more dissipative since the pseudo-thermal buoyancy force increases in 
the same limit. Hence, the RBNP is unlikely to occur in the tracer limit for moderately 
dissipative particles. 

To probe the effect of the volume fraction of the intruders, we show the phase-diagram for 
the BNP/RBNP transition in Fig. [6]for three relative volume fractions: <fii/<f)s — 10~ 8 , 0.1 
and 1. The other parameters are as in Fig. [2J As we increase the relative volume fraction of 
the intruders, the ranges of mass- and size-ratios for which the reverse segregation (RBNP) 
occurs increase sharply. It is interesting to recall that most of the experiments of Breu 
et al. [kJ correspond to the case where the number of layers of each species was equal; 
this translates into <pi > <p s . However, when they reduced the number density of larger 
particles they could not observe RBNP (with other conditions remaining the same). Also, 



the recent experiments of Yan et al. 



371 ] (at reduced air-pressures) do not show RBNP. Our 



model predictions are, therefore, in qualitative agreement with these experimental findings. 
Clearly, more experiments are needed to map out the correct phase-diagram for various 
volume fractions of the intruders. 



V. DYNAMICS OF INTRUDERS IN TRACER LIMIT: NONMONOTONIC 

RISE TIME 

As an application of our theoretical framework, we now consider the dynamics of intruders 
in the tracer limit (ni << n s ) and determine the intruder's relative velocity u\ and its rise- 
time. (Hereafter, the dynamics of intruder particles in the tracer limit is considered to be 
that of a single intruder in a bed of small particles.) Note that these two quantities have 

nnn 

been measured [5|, [6|, 139|] in many Brazil-nut experiments. 

To probe the motion of the intruder, we need to know the macroscopic velocity-field u s 
of the smaller particles a priori. Let the system be excited by a periodic force in the vertical 
direction with a (symmetric) harmonic displacement: z p (t) = Asm(ut). We assume that 
the vibrofluidized state of the smaller particles is coupled to this periodic movement, and 
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make the following approximation for the velocity u s of the smaller particles: 

dz 

u s « = -Ajcos(cut). (43) 

This simply implies that the smaller particles follow the motion of the base plate which is a 
reasonable assumption for shakings at a high frequency but with low amplitudes (T >> 1). 

We need to estimate the granular temperature of the mixture T which appears in the drag 
term in the evolution equation ( |23i) . In the tracer limit (n; << n s ), the mixture granular 



temperature would be that of an equivalent mono-disperse system of smaller particles 33]; 
we have already assumed that the granular temperature is uniform throughout the bed. 
Considering this homogeneous-state of ffuidization, the global temperature can be estimated 
by equating the rate of energy input through the bottom plate with the rate of energy loss 
due to inelastic particle collisions [38] : 



1 m s S p {Auf 
V2n Nr* (1 - ' 1 ' 

Here S p is the surface area of the base-plate and N is the total number of particles. As 
demonstrated in many earlier studies, this provides a reasonable approximation for the 
avera ge t emperature of the bed, and represents the leading-order solution even in the dense 
limit [38]. We assume that in the fluidized state the height of the bed is H. Thus, the 
expression for the global temperature can be rewritten as 

= V2 m s A 2 u 2 

3 4>sRhRis{1 - el)' 1 ' 

where 

R " = J; W 

is the non-dimensional bed-height as a multiple of the intruder-diameter, and Ri s = r\jr 8 
the size-ratio. 

Let us introduce the following reference scales for non-dimensionalization: 
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CO' 1 



(47) 
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where the quantities with overbars are non-dimensional. Hereafter, for convenience, we will 
drop the overbar on the non-dimensional quantities. With this scaling, the non-dimensional 
granular temperature has the following expression: 

rp _ (48) 

3<I>.RhRi.(1 - e 2 „Y K } 

The evolution equation ff23l) in non-dimensional form can now be written as 

-jp- = a (r _1 — sint) — (3u\ (49) 

where 

(3 = (3 R A VT (50) 

is the non-dimensional drag coefficient which is a function of the bed-temperature T and 
the amplitude of the harmonic-shaking 

= < 51 > 

Here a and (3q are non-dimensional functions of the size-ratio, mass-ratio and volume frac- 
tion: 

/ m s Z t Ti _ ^ 
V m Z s T s 

a = sk„ ( ( i^ M : 1/2 

\Ris + I J \ir mi 

K U = ^sXlsO- + Rlsf 
m s 

M sl = 

mi + m s 
r i 

Ris = —■ 

r s 

The solution to the differential equation (1491) . with initial condition u[{t)\ t= Q = 0, is given 

by 

u\{t) = jjL (1 - C -*) - (jfp) - cost + (3 sin t) . (52) 
The temporal evolution of the position of the intruder, zi(t), can be obtained by integrating 
Eq. (J52D: 

with zq = ^(0) being the position of the intruder at time t — 0. Note that has been 
non-dimensionalized by the diameter of the intruder. 
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A. Rise Time: Non-monotonicity and Experiments 



By knowing the steady relative velocity of the intruder w[, we can calculate its rise-time 

T = ^i = ^K. (54 ) 
u\ a 

This is the (non-dimensional) time that the intruder will take to travel a height of one 
intruder diameter. In the following we have plotted the rise-time in cycles. 

First, we consider an intruder of a given size-ratio (Ri s = r\jr s = 2.5) and vary its density; 
further, we assume the equipartition assumption TJ = T s . For this case, the variation of the 
rise-time with the density-ratio, R p = pi/p s , is shown in Fig. [7b; other parameters are as 
in Fig. HI The parameters for base-plate excitations are T = 5 and / = 50 Hz; we set the 
bed-height to Rh = H/{2r{) = 10. As expected, r diverges at the onset of BNP/RBNP 
(R p i ~ 0.64) since the intruder velocity, u[ changes its sign at this threshold density as seen 
in the inset of Fig. [7b. For R p > R p \, the sink-time decreases, implying that the heavier 
intruders sink faster; for R p < R pl , the lighter intruders rise faster. 

Now we relax the equipartition assumption, i.e. T} ^ T s , and plot the variation of r 
with R p in Fig. [7b In contrast to the case in Fig. [7b, now we have the pseudo-thermal 
buoyancy force acting on the intruder. This additional upward force makes the intruder 
to rise again beyond R p > R p2 . Interestingly, the rise-time for R p > R p2 decreases with 
increasing density-ratio, implying that the heavier the intruder the faster it rises. Such an 



effect has been observed in experiments 12J, |37J, [39J. Within the window of the RBNP 



(i? p i < R p > Rpz), we observe that the sink-time varies non-monotonically and diverges at 
both the limits. For this parameter combination, an intruder with R p = 1 will take about 
15 cycles of excitation to travel (sink) a distance of its diameter. 

With other parameter conditions as in Fig. ©, we increase the intruder size such that we 
are in the zone of BNP (see Fig. [1]). This is shown in Fig. [8] for three different size-ratios 
(Ris = 3,4,5). It is interesting to note that the rise-time is non-monotoic with the density- 
ratio. For a given size-ratio, there is a threshold density-ratio (R p < 1) for the intruder 
above/below which it rises faster. This threshold density-ratio shifts to a lower value with 
increasing size-ratio. Note further that at a given density-ratio the larger the intruder the 
faster it rises in conformity with our earlier observation in Fig. [7f>. 

Next we fix a large size-ratio of Ri s = 10, and show the variation of r with the density ratio 
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R p in Fig. [9] for four different restitution coefficients (e = 0.999,0.998,0.997,0.99). With 
increasing dissipation, the magnitude of the energy-ratio, T}/T s , increases that generates 
additional pseudo-thermal buoyancy force on the intruder, thereby increasing its velocity. 
Hence the intruder will move faster and the height of the peak in the rise-rate curve decreases 
as observed in Fig. [9J 



B. Comparison with Experiments 

In order to verify the possible effects of convection J4], here we make a qualitative com- 
parison of our predictions on the rise time with the recent experiment results of Huerta 
and Ruiz-Suarez 39j. They performed two-sets of experiments to test the effects of bulk 



convection, one at / = 5Hz and T = 3 (i.e. the low frequency but high amplitude limit), 
and the other at f — 50Hz and T = 3 (i.e. the high frequency but low amplitude limit). 
(They verified that the bulk-convection was negligible in the high frequency limit, but was 
present in the low frequency limit.) They observed that the rise-time is non-monotonic 
in the high-amplitude (and low frequency) limit (see their Fig. 1), but in the absence of 
convection the heavier particles sank (see their Fig. 3). 

Note that most of the assumptions in our model are in tune with their high frequency 
experiments since the convection effects are claimed to be negligible in that limit. However, 
if the level of pseudo-thermal buoyancy force is small, the window of the RBNP (R p i < 
R p < R P 2, refer to Fig. [7t>) can be made arbitrarily large as in Fig. 02. In this case, the 
Archimedean buoyancy force competes with geometric forces, thereby deciding the dynamics 
of the intruder for R p > R p \ as in our Fig. [Tfct which looks remarkably similar to Fig. 3 in 



ref. 



39]. 



Moving onto their low frequency experiments, we need to consider the additional effects of 
bulk-convection. Note that the bulk-convection generates an upward-plume at the center of 
the container, and two 'narrow' downward plumes near the side-walls; clearly, the granular 
energy of the particles within the center-plume will be much greater than that of bulk 
material, and this would lead to another competing force on the intruder. Under the present 
formalism, this additional force would reinforce the pseudo-thermal buoyancy force. Thus, 
the effective pseudo-thermal buoyancy force will increase in the presence of convection, 
leading to an over-turning of the phase-coexistence line beyond a critical size-ratio (refer to 
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our Fig. [Q. Hence the added effect of bulk-convection is to take us to the regime of BN 



and the heavier particles will also rise. This is precisely what we see in Fig. 1 of ref. [39]. 
Therefore, in their low frequency experiments, the effective pseudo-thermal buoyancy force 
would compete with the Archimedean buoyancy force and geometric forces. 



Lastly, we comment on the experiments of Mobius et al. [12] who first showed that 
the rise-time of an intruder (rj/r a > 10) in a bed of small particles (r s = 250/xm) is non- 
monotonic with the density-ratio. They also showed that the peak on the rise-time curve 
decreases in height and shifts to a lower density with decreasing the air-pressure, and vanishes 
as the air-pressure approaches 1 torr. It is conceivable that the decreased air-pressure has 
reduced the effective drag on the intruder and thereby increased its velocity- hence the 
height of the peak (maximum rise-time) decreases with decreasing air-pressure. But the 
shifting of this peak to a lower density remains unexplained. 

VI. DISCUSSION 

A. A Simple Explanation for RBNP: Boltzmann Limit 

The reverse segregation effect can be understood by considering the Boltzmann-limit (i.e. 
the dilute limit </> ^ 0) for which the equation of state of species i is 

Pi = riiTi. 

By integrating the steady vertical momentum equation, we obtain an expression for the 
number density profile of species i, 



m{z) = m(0) exp 



rriigz 



(55) 



T 

where we have assumed that Tj is independent of z. (This assumption implies that the 
variation of number- density with z is due to the variation of partial pressure along z.) 
Hence, the ratio of number densities is 



m(z) rii(0) 

exp 



Ti V T s mi 



(56) 



n s {z) n s (0) 

For a mixture of particles with purely elastic collisions, the equipartition principle holds, 
i.e., Ti = T s . In this case, the particles segregate according to their masses, driven by the 
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Archimedean buoyancy force, and there is no segregation if the particles are of the same 
mass. Thus, the heavier particles remain at the bottom and lighter particles at the top, 
leading to the RBN effect. 

Now considering inelastic particles of equal masses m; = m s but having different sizes 
and densities, we know that TJ < T s for all values of restitution coefficient. This implies 
that the decay-rate of n\ with z is faster than that of n s . Hence the center of mass of the 
larger-particles will be at a lower-level than that of the smaller-particles, 

^ %l ^ 'cm -, 
~~ ~ - ' , 

^ %s cm 

clearly leading to the RBN effect. 



B. Higher-order Effects and Possible Experiments 

The granular hydrodynamic model is a natural and convenient framework for studyin 



fluidized systems. It is nonetheless criticized for being overly simplified and unrealistic [14]. 
In our analysis we have considered the hypothesis of a Maxwellian velocity distribution 
function. Let us remark that the leading-order solution of the Boltzmann-Enskog kinetic 
equation for granular mixtures is a Maxwellian (see ref. [201]), and the non-Gaussian cor- 
rection term remains relatively small ]31]. More importantly, since the present theory is 
restricted to the Euler-level description of hydrodynamic equations, we do not need to take 
into account the corrections due to non-Maxwellian effects as explained in jl9]. Moreover, 
the Archimedean buoyancy force ) and the static geometric force {F^) do not depend 
on the distribution function, and the non-Gaussian correction does not affect the dynamic 
geometric force (F^ n ). Only the pseudo-thermal buoyancy force (-F 6 T ) can be expressed in 
terms of the velocity distribution function. Thus, the higher-order corrections would affect 
this buoyancy force- a detailed investigation of this is left for a future investigation. 

All the segregation forces are indirectly measurable in a standard vibrated-bed setup, 
and thereby making them directly verifiable via experiments. For example, the thermal 
buoyancy force {F^f) can be measured by measuring the granular energies (7] and T s ) from 
the snapshots of successive particle configurations. These snapshots can also be used to 
measure the pair correlation function (xij an d hence the compressibility factor Z,) and the 
dynamic geometric force (F^ n ). The other two forces (F^ and i* 1 **) are then easily obtained. 
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Knowing all these forces from experiments, one could construct each term of our theory 
independently. We hope that the present investigation will stimulate future experiments to 
measure these quantities. 

VII. SUMMARY AND CONCLUSIONS 

One of our motivations has been to explore whether the hydrodynamic-description of 
granular materials can be applied to the derivation of a time-evolution equation for the 
segregation velocity of intruders in a dense fluidized bed. We have done that in our previous 
paper [3] , starting from the kinetic-theory of binary mixtures. Based on this time-evolution 
equation, here we have tried to provide a unified theoretical description for the segregation 
dynamics of Brazil-nuts in a vibrated granular mixture. 

We have discussed a novel mechanism for segregation and argued that the onset of segre- 
gation is due to a competition between the buoyancy and geometric forces. Apart from the 
standard Archimedean buoyancy force, we have introduced the notion of a pseudo-thermal 
buoyancy force that results from the fact that the intruder particle fluctuates at a different 
energy level than the smaller particles. (The microscopic dissipation, due to the inelastic na- 
ture of particle-collisions, is responsible for the non-equipartition of granular energy.) The 
size-difference between the intruder and the bed-material results in two geometric forces, 
one compressive and the other-one tensile in nature- the net geometric force is always com- 
pressive. Thus, the net buoyancy force competes with the net compressive geometric force, 
thereby deciding the onset of BNP/RBNP transition. 

For a mixture of perfectly hard-particles with elastic collisions, the pseudo-thermal buoy- 
ancy force is zero but the intruder has to overcome the net compressive geometric force to 
rise. For this case, the competition between the net geometric force and the Archimedean 
buoyancy force yields a threshold density-ratio, R p \ = pi/p s < 1, above which the lighter 
intruder sinks, thereby signalling the onset of the reverse buoyancy effect. For a mixture of 
dissipative particles, on the other hand, the non-zero pseudo-thermal buoyancy force gives 
rise to another threshold density-ratio, R P 2 ( > R p i), above which the intruder rises again. 

This theory is applied to study the dynamics of a single intruder (i.e. in the tracer limit 
of intruders) in a vibrofluidized granular mixture, with the effect of the base-plate excitation 
being taken into account through a 'mean-field' assumption. The rise (/sink) time diverges 
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near the threshold density-ratio for reverse-segregation. The most interesting result is that 
the rise-time of the intruder could vary non-monotonically with the density-ratio. For a 
given size-ratio (in the zone of BNP), there is a threshold density-ratio for the intruder at 
which it takes maximum time to rise and above(/below) which it rises faster. This implies 
that the heavier (and larger) the intruder, the faster it ascends; a similar effect has been 



observed in some recent experiments 12|, |39|. The peak on the rise-time curve decreases in 



height and shifts to a lower density-ratio as we increase the magnitude of the pseudo-thermal 
buoyancy force. 

The main message emerging from our analysis is that the pseudo-thermal buoyancy force 
that results from the non-equipartition of granular energy (T} ^ T s ) plays an important 
role in the segregation process. Even though the source of this energy non-equipartition is 
inelastic dissipation, it can be argued that the effects of the bulk-convection will also lead to 
a separation between the granular energies of the intruder and the bed-particles even in the 
limit of perfectly elastic collisions. Thus, the increased pseudo-thermal buoyancy force (i.e. 
with increasing dissipation-levels in our model) can also be related to the presence of the 
bulk-convective motion. In this connection it may be pointed out that a different functional 
form for the equation of state might change the precise threshold of segregation, but the 
presence of additional pseudo-thermal buoyancy force, for example, due to convection, is 
likely to retain the overall phase-diagram similar. 

It is remarkable that our theory can explain several experimental observations in a unified 
manner, despite having many simplifying assumptions. In addition to explaining the exper- 
imental results of Breu et al. [3| on the reverse Brazil-nut segregation, we have suggested 
a plausible explanation for the origin of the onset of the reverse-buoyancy effect of Shinbrot 
and Muzzio [l0(] as well as the origin of the non-monotonic ascension-dynamics 12j, |39 |. 
However, the effects of air-pressure on the non-monotonic ascension- dynamics [l2j seems 
to be more subtle and needs further work. This can be done by considering a three-phase 
mixture, with separate balance equations for the air and taking into account the resulting 
interactions on both the smaller particles and the intruders. Lastly, the effects of Coulomb 
friction and the particle roughness need to be incorporated in our model in future. 
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FIG. 1: (color online) Effect of inelasticity on the phase diagram for BNP/RBNP in three dimen- 
sions: <f) = 0.5 and <t>i/(j>s = 10~ 8 . The dotted line represents a mixture of equal density particles 
(Pi = Ps)- mi/m s = (n/r s ) 3 . 
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FIG. 2: (color online) Variations of segregation forces (F/m s g) with the size-ratio for a mixture of 
equal-density particles (pi = p s ) at e = 0.95. The dotted arrow indicates the locus of the transition 
BNP 44> RBNP. The upper inset shows the variations of the static and dynamic contributions to 



the total geometric force {F t 
of Ti/T s with the size-ratio 



dyn^ 



26]. 



TP st _i_ TP' 



with the size-ratio. The lower inset shows the variation 
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FIG. 3: (color online) (a) Variations of segregation forces (F/m s g) with the size-ratio for a 

mixture with fixed mass-ratio (mi/m s = 10) at e = 0.95. The dotted arrow indicates the locus 

of the transition BNP 44> RBNP. The upper inset shows the variations of the static and dynamic 

contributions to the total geometric force (F ge = F^ e + Fge n ) with the size-ratio. The lower inset 

shows the variation of Ti/T s with the size-ratio, (b) Variations of Archimidean (F^) and thermal 

39 

(Fjf) buoyancy forces with ri/r s , with parameter values as in (a). The arrow on the i*j -curve 

;„j;„„j-„„ j-l,„ „;„„ „„-!-;„ rpT ^ n 




FIG. 4: (color online) (a) Onset of reverse buoyancy effect in terms of segregation forces: variations 
of F/m s g with the mass-ratio for a mixture with ri/r s = 2.5. Parameter values are <ft = 0.5, 
4>i/4> s = 10 -8 and e = 0.95. The vertical arrows indicate the loci of the transition BNP 44> RBNP. 
The upper inset shows the variations of the Archimedean and pseudo-thermal contributions to the 
total buoyancy force (F^ = + F?) with the mass-ratio. The lower inset shows the variation of 
Ti/T s with the mass-ratio, (b) Same as in a but ^ith density-ratio, pi/p s - The redline in the lower 
inset shows the variation of F^ + F ge . 
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FIG. 5: (color online) Effect of mean volume fraction on the phase diagrams for BNP/RBNP 
in three dimensions: e = 0.95 and 4>ll4>s = 10~ 8 . The dashed arrow indicates the direction of 
increasing shaking strength T. 
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FIG. 6: (color online) Effect of the relative volume fraction of intruders on the phase diagram for 
BNP/RBNP in three dimensions: eft = 0.5 and e = 0.95. 
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FIG. 7: (color online) (a) Variations of rise (/sink) time r (in cycles) with density-ratio for T\ = T s ; 
the parameter values are as in Fig. |H The shaking parameters are T = 5, / = 50 Hz; the intruder 
diameter is 4.97cm and the bed-height Rh = H/{2r{) = 10. (b) Same as a with non-equipartition 
assumption. The inset in each panel shows the variation of the intruder's relative velocity, vX, 
with the density-ratio. The vertical lines in eachj^anel refer to density-ratios for the BNP/RBNP 
transition. 




FIG. 8: Variation of rise time r (in cycles) with density-ratio for different size-ratios; the parameter 
values are as in Fig. [5j The shaking parameters are as in Fig. 
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FIG. 9: Variation of rise time r (in cycles) with density-ratio for different restitution coefficients. 
The size-ratio is set to ri/r s = 10; other parameters as in Fig. 
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